Effect of anharmonicity of the strain energy on band offsets in semiconductor nanostructures 
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Anharmonicity of the inter-atomic potential is taken into account for the quantitative simulation of the con- 
duction and valence band offsets for highly-strained semiconductor heterostructures. The anharmonicity leads 
to a weaker compressive hydrostatic strain than that obtained with the commonly used quasi-harmonic ap- 
proximation of the Keating model. Inclusion of the anharmonicity in the simulation of strained InAs/GaAs 
nanostructures results in an improvement of the electron band offset computed on an atomistic level by up to 
100 meV compared to experiment. 
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The accurate simulation of the electronic structure is of 
utmost importance for the design of nanoelectronic and op- 
toelectronic device structures. It has been shown both 
theoreticallyit and experimentalljiliili^, that the energy spec- 
trum in semiconductor nanostructures is extremely sensitive 
to the built-in strain. Pryor et aim has shown that the contin- 
uum elasticity method fails to adequately describe the strain 
profile in InAs/GaAs heterostructures with a 7% lattice mis- 
match between the constituent materials. The two-parameter 
valence-force-field (VFF) Keating model 6 - 7 is a commonly 
used approximation for atomistic-level calculations of the 
equilibrium atomic positions in realistic-size nanostructures 8 . 
This paper shows that the quasi-harmonic Keating model is 
insufficient to describe highly strained InAs/GaAs nanostruc- 
tures, where anharmonicity of the strain energy becomes im- 
portant. 

Keating's model treats atoms as points connected with 
springs in a crystal lattice. The strain energy depends only 
the nearest-neighbor interactions within a quasi-harmonic 
approximation 
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The first coefficient, a, corresponds to the spring constant for 
the bond length distortion, while the second one, j3, corre- 
sponds to the change of the angle between the bonds or, so- 
called, "bond-bending". Here the summation is over all m 
atoms of the crystal and their nearest neighbors n and k. r mn 
and d mn are the vectors connecting the m-th atom with its 
n-th neighbor in the strained and unstrained material, respec- 
tively. 

A realistic interatomic potential is schematically shown in 
Fig. ^ with a line marked by large circles. The Keating in- 
teratomic potential corresponding to the Eq. Q is plotted in 
Fig. \l\ with dashed lines. The Keating potential is referred 
to as "quasi-harmonic" rather than "harmonic" since it de- 
pends upon the fourth rather than the second power of the 
bond length. This quasi-harmonic potential fails to reproduce 
the weakening of the interatomic interaction with increasing 
distance between atoms and it underestimates the repulsive 
forces at close atomic separation. Therefore Eq. Q can ad- 
equately describe the strain energy only at small deforma- 
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FIG. 1: Schematic interatomic potential used in the Keating (dashed 
line) and our model (solid line). Dash-dot line plot the model in- 
teratomic potential with the anharmonicity corrections to the VFF 
constants before the truncation. The line marked with large circles 
approximafly traces the shape of the realistic interatomic potential. 



tions. In the widely used InAs/GaAs heterostructures, the 
lattice mismatch is as large as 7% and anharmonicity of the 
interatomic potential is expected to become important. 

We include the anharmonicity of the interatomic poten- 
tial directly into the Keating model by introducing distance- 
dependent VFF constants a and j3: 
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B, 

,C rf . Q nmk and Ql mk are 
the actual and the unstrained angles between mn and mk 
bonds, respectively. In homogeneous mono-atomic and bi- 
nary compounds, all bonds are the same and the indexes to, 
n, and k can be dropped. a and (3 are the VFF constants in 
the unstrained material. The anharmonicity correction coef- 
ficients A and C describe the dependence of bond-stretching 
and bond-bending forces on hydrostatic strain, while B is re- 
sponsible for the change of the bond-bending term with the 
angle between bonds. The details of the derivation of the 
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anharmonicity corrections A, B, and C from the experimen- 
tal phonon spectra of strained bulk materials are presented in 
Ref A The parameters used for the simulation of InAs/GaAs 
nanostructures are listed in Table|I] 



Material 


a (N/m) 


A) (N/m) 


A 


B 


C 


GaAs 


41.49 


8.94 


7.2 


7.62 


6.4 


InAs 


35.18 


5.49 


7.61 


4.78 


6.45 



TABLE I: Valence-force-field constants in unstrained materials and 
anharmonicity correction coefficients for InAs and GaAs. 

Introduction of the anharmonicity corrections in the VFF 
expression makes the form of the model potential more real- 
istic and expands the validity range of the strain simulations 
(see Fig. [Q. In order to ensure the convergence of the min- 
imization of the energy Q with a and (3 given by Eqs. (|2ji, 
l|3}, our model inter-atomic potential was truncated as shown 
in Fig. ^ with the solid line. Note that the proposed introduc- 
tion of the anharmonicity corrections in VFF constants does 
not substantially increase the computational time of the mini- 
mization of the total strain energy, since it does not introduce 
any additional summation in Eq. Q. 

To illustrate the effect of the anharmonicity on the strain 
distribution in III-V semiconductor nanostructures, the hy- 
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l/6( £ xa: + £yy — ^zz), components of the strain in ultra-thin 
epitaxial InAs/GaAs layers have been computed using both 
the conventional Keating model and our model. The geom- 
etry of the structures and the obtained built-in strain are pre- 
sented in Table lU Comparing the results of the two mod- 
els, we note that the sharp rise of the strain energy at small 
inter-atomic distances leads to a smaller equilibrium hydro- 
static compression than is obtained with the Keating model. 
At the same time, the quasi-harmonic approximation underes- 
timates the bond stretching deformation. In contrast, biaxial 
compression is increased in our anharmonic model, while the 
biaxial tension is suppressed. 

The band offsets for different InAs/GaAs nanostructures 
obtained for the strain distribution simulated within the Keat- 
ing's model and including the anharmonicity corrections are 
compared with the available experimental data 3,4 ' 5 in Tablelllll 
The local band structure was obtained within the sp 3 d 5 s* em- 
pirical tight-binding model where the Hamiltonian matrix el- 
ements depend on the distance between the atoms^l The 
tight-binding parameters were fitted to reproduce the proper- 
ties both of the unstrained and strained bulk materials. Ultra- 
thin epitaxial layer structures^ provide us with the perfect 
opportunity to test the effect of the anharmonicity at strong de- 
formation on the valence and conduction band offsets. Anhar- 
monicity corrections reduce the discrepancy between the ex- 
perimental and modeled energies significantly (see Table IIIII. 

The strain distribution (Fig. [2] a, b) and the energy spectrum 
(Fig. |2 c) were computed for the quantum dot crystal (QDC) 
reported in Refm. The structure consists of three layers of reg- 
imented dome-shaped quantum dot arrays (with a 20 nm base 
diameter and a 7 nm height), vertically stacked with a small 
(3 nm) vertical separation between the QD layers. The 0.7 nm 
wetting layer is also included in our simulations to properly 
model the electronic spectrum' 1 . The built-in strain distribu- 



CB 
i 



n5 
o 



in 
o 
i— 

x ■ 



r 






(A 


r 


InA 




InA 


iAs 




lAs 





Wetting 
layer 



5 
H 

X 

s 
in 



cn 
i— 

CD 
c 
CD 



£ 
o 



o 

JD 
LU 




20 40 60 

Position along the growth direction (nm) 

FIG. 2: Computed distribution of the hydrostatic (a) and biaxial (b) 
strain components, and (c) electronic band structure along the growth 
direction in InAs/GaAs QDC. The structure taken from Ref«« consists 
of three layers of dome-shaped QDs with a 20 nm base and a 7 nm 
height on top of the 0.7 nm wetting layer. The crossection is made 
near the center of the quantum dot stack. The results obtained with 
the Keating model are plotted with black dots. The results obtained 
with the anharmonic model are plotted with solid line on (a) and (b) 
and with gray dots on (c). The thin lines on (c) show the edges of 
conduction, valence, and spin-orbit split-off bands at the center of 
the Brillouin zone in the unstrained materials. 



tion in such structures is very inhomogeneous. The average 
hydrostatical component of the strain tensor inside InAs quan- 
tum dots (Fig. |3 a) is overestimated by about 25% within the 
commonly used Keating model. At the same time, the biax- 
ial strain distribution (Fig. [2 b) changes little when computed 
with the different models. 

The main effect of the anharmonicity is a downward shift of 
the conduction band edge inside the quantum dots, as shown 
in Fig. |2 c. This is caused by the sensitivity of the conduction 
band to the hydrostatic compression, which is smaller within 
the anharmonicity model. The difference in the correspond- 
ing band edges obtained within the two models is as large 
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Ref. 




Structure 


Size 


Substr. 


Hydrostatic strain (%) 


Biaxial strain (%) 








L x 


Ly 






4 


S H 




A r 

e B Ob 




1.0 0.0 


SQW 


2 ML 


5 ML 


GaAs 


-2.97 


-2.59 


14.9 


-3.71 


-4.10 -9.4 




0.0 1.0 


SQB 


2 ML 


5 ML 


InAs 


3.54 


4.30 


-17.7 


3.72 


2.90 28.3 




1.0 0.0 


MQW 


1 ML 


30 nm 


GaAs 


-2.87 


-2.33 


23.5 


-3.81 


-4.38 -13.0 



TABLE II: Hydrostatic (H) and biaxial (B) strain components within the thin epitaxial layer in several structures computed within the Keating 
(K) model and taking into account the anharmonicity of the strain energy (A) for different In^Gai-ajAs/InyGai-yAs nanostructures^-. Nota- 
tions: the nanostructure type: SQW - single quantum well, SQB - single quantum barrier, MQW - multiple quantum wells, ML - monolayer, 
L x - width of the In^Gai-^As epitaxial layer, L y - width of the LiyGai-yAs epitaxial layer (corresponds to the thickness of the capping 
layer for SQW and SQB). 5 — (e K — e A )/e A shows the relative change of the computed built-in strain when the anharmonicity is neglected. 



Ref. 


Structure 


Experimental 
method 


K 


Sk (%) 


AE C (meV) 
A 5 A (%) 


Exp. 


K 


5k (%) 


AE V (meV) 
A 5 A (%) 


Exp. Sub-band 




SQW 


XPS 


471.5 




574.0 




374.6 


-34.4 


447.0 -15.7 


530 ± 50 hh 




SQB 


XPS 


-40.7 




-91.3 




-271.7 


69.8 


-225.9 41.2 


-160 ±50 Ih 




MQW 


CV&DLTS 


475.7 


-31.1 


584.4 -15.3 


690 


373.4 




430.5 


hh 




QDC 


DLTS 


242.0 


-29.0 


347.0 1.8 


341 ± 30 











TABLE III: Experimental band offsets in the conduction (AE C ) and valence (AE V ) bands compared with the offsets computed within the 
sp 3 d 5 s* empirical tight-binding model— using the equilibrium atomic positions found within the two-parameter Keating model (K) and 
including anharmonicity corrections to the VFF constants (A) for different In^ Gai As/In,, Gai_ y As nanostructures. The SQW, SQB, and 
MQW structures are the same as in Table ItTI The QDC structural consists of three layers of dome-shaped QDs with a 20 nm base and a 7 nm 
height on top of a 0.7 nm wetting layer. The distance between the QD layers is about 10 nm. The band offsets are determined as the following: 
AE C = i? c (In y Gai_ y As) — ^(In^Gai-^As) and AE V = — (^(InyGai-^As) — ^(In^Gai-a; As)) so they would be positive for potential 
well and negative for potential barrier. Notations: the nanostructure type: SQW - single quantum well, SQB - single quantum barrier, MQW 
- multiple quantum wells, QDC - quantum dot crystal, i.e., the 3-dimensional ordered array of quantum dots; the band offset measurement 
techniques: XPS - X-ray photoemission spectroscopy, CV - capacitance-voltage spectroscopy, DLTS - deep-level transient spectroscopy. 5k 
and 5a estimate the relative deviations of the simulation from the experiment. 



as 105.0 meV. This shift brings the overall band offset be- 
tween the InAs quantum dot and GaAs buffer computed in our 
model very close to the experimentally observed value (see 
Table flTO. Due to the small difference in the biaxial strain 
distribution obtained with the two models (see Fig. [2] b), the 
energy structure of the valence band remains almost the same, 
as can be seen in Fig. [2 c. 

In conclusion, we have demonstrated that the anharmonic- 
ity is important for the modeling of the electronic states in 
highly strained InAs/GaAs system. Compared to the standard 
Keating model we have found corrections of over 100 meV 
in some band offsets, resulting in values significantly closer 
to the experimental data. This demonstrates that the deforma- 



tion in the nanostructures is beyond the applicability range of 
the quasi-harmonic approximation for the strain energy. The 
anharmonicity corrections can be performed without a signif- 
icant increase of the computational cost, since the model re- 
mains limited to a nearest neighbour interactions. 
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